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o 
o 

, We consider the analytical properties of the single soliton solution in a Skyrmion-type 



Oh 



X 



Lagrangian that incorporates the scaling properties of quantum chromodynamics (QCD) 



■ through the coupling of the chiral field to a scalar field interpreted as a bound state of 
gluons. The model was proposed in previous works to describe the Goldstone pions in a dense 
medium, being also useful for studying the properties of nuclear matter and the in-medium 
properties of mesons and nucleons. Guided by an asymptotic analysis of the Euler-Lagrange 

■ equations, we propose approximate analytical representations for the single soliton solution 

in 
o 
^i- 
o 

pH , We find that the convergence of the approximate representations to the numerical solutions 

is considerably improved by taking the expansion coefficients as free parameters and then 



in terms of rational approximants exponentially localized. Following the Pade method, we 
construct a sequence of approximants from the exact power series solutions near the origin. 



minimizing the mass of the Skyrmion using our ansatze for the fields. We also perform an 
analysis of convergence by computation of physical quantities showing that the proposed 
analytical representations are very useful for further phenomenological calculations. 

I. INTRODUCTION 

Through the formalism of effective field theory a significant progress has been reached in our attempt 
to undertstand the physics involving mesons and baryons in nuclear matter. The first effective descrip- 
tions appeared as generalizations of the linear a model so that the physical constraint dictated by the 
broken-realized chiral symmetry of QCD was naturally included. More recently, the interest was focused 
on models incorporating the broken scale invariance of strong interactions and therefore suitable for 
describing physics beyond the pion sector Q, Qj. Among these models, an interesting approach consists 
in introducing a scalar glueball with an appropriate potential into a nonlinear realization of the a model, 
thus accounting for both broken-scale and chiral symmetries From this background, the authors in 
Ref. jj| proposed a Skyrmion-type Lagrangian with spontaneously broken chiral and scale symmetries. 
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This model, which is valid as an approximation to QCD in the large iV c -limit, is useful for studying the 
properties of pions in dense medium and can also describe multibaryon systems, including dense hadronic 
matter. In this scenario, baryons arise as solitonic solutions of the equations of motion that, because of 
their nonlinearity, cannot be solved in an analytical closed form. 

Naturally, it is desirable that the advance in the theoretical domain be accompanied by the solution of 
the mathematical problems introduced. Our contribution in this paper is to give analytical approximate 
representations for the single soliton that is used to define the parameters of the theory at zero density. 
This types of representation are preferred to numerical solutions because with them one can trace exact 
analytical information about the soliton solution. Besides, explicit representations are very convenient 
for phenomenological calculations in the analysis of dense Skyrmion matter. Our procedure will consists 
of extracting the structure of the exact solution from an asymptotic analysis of the equations of motion 
for the scalar and chiral fields, and then we look for solutions in terms of rational approximants, Pade 
Approximants, that have proved to be particularly successful in giving a faithful representation of the 
topological configuration of the original Skyrme Model [5|. Here, the task is more involved since it 
requires the resolution of two coupled nonlinear differential equations for the new scalar field and the 
chiral field already present in the Skyrme Model. Furthermore, the nonzero mass for pions and the 
additional scalar field imply exponential asymptotic behaviors for both solutions, and this feature must 
be incorporated in the approximants. We will see that the analytical approach to the single Skyrmion 
with the scalar field problem can be implemented in a simple way and that the suggested representations 
reproduce the properties of the exact solutions for small and large values of the radius r. 

Our paper is organized as follows. In the next section we briefly describe the model that is composed 
by the original Skyrme Model implemented with the trace anomaly of QCD, signal of the broken-realized 
scale invariance. In Sec. Ill we present the properties of the single Skyrmion and the scalar field and then, 
in Sec. IV we construct the rational representations by using the exact asymptotic solutions for both the 
chiral and the scalar field and we present an analysis of convergence of the proposed approximants. 
Finally, we give the conclusions of our work in Sec.V. 



II. THE MODEL 



The model constitutes an extension to the original Skyrme Model that emerged as an attempt to 
describe baryons interacting by means of meson exchanges 6]. In a more recent formulation introduced 
by Ellis and Lanik [3] and developed in further works for the description of nuclear matter and finite 
nuclei [||, the task was aimed at incorporating the underlying scale invariance of QCD under the scale 
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transformation 

x->\~ l x, A>0. (1) 

This symmetry is broken at the quantum level by dimensional transmutation, and its breaking is 
manifested through a nonvanishing trace of the energy momentum tensor 9^, referred to as the QCD 
trace anomaly which in the chiral limit is given by 

0£ = ^TrG^CT, (2) 

where /3(g) = -(g 3 /48vr 2 )(lliV c - 27V/) is the one-loop beta function of QCD. 

The description starts with the Skyrme Lagrangian for massive pion fields, which reads 

£ = llTr(d^d^U) + J-Tr[C/t«9 - U, U^8 U U} 2 + &^Tr(U + - 2), (3) 
4 32e z 4 

where the Goldstone pion fields 7r a play the role of a chiral phase angle present in the chiral field 
U = exp(ir • tt/ f n ) belonging to SU(2), and /„-, and e are the pion decay constant, the pion mass 
and the Skyrme parameter, respectively. 

The broken-scale invariance of QCD is implemented by the introduction of a scalar field x( r ) with a 
scale dimension 1. This scalar degree of freedom is interpreted as a bound state of gluons, the glueball 
field, that operates on two-pion exchange physics in meson-exchange phenomenology. 

With an appropriate scaling breaking-potential for the scalar field x{ r )i the modified Lagrangian that 
respects the scaling properties of QCD reads 



C = & (^j Tr(d^d»U) + J-^Tr([C/t^[/, U^d v U]f + f -^- (j^j Tr(U + rf-2) 

(4) 



1„ „„ lm? 



Here, f x is the vacuum expectation value of the scalar field x( r )- The mass of x( r ) 1S given by m 2 = 
d 2 V(x)/dx 2 , where V(x) denotes the scale-breaking potential. No experimental values are available for 
those parameters, and conjectural values depend rather on specific considerations of the model, as it is 
explained in Ref. Q]. The inability to fix them comes from the lack of knowledge about the mixing 
between the two components of the trace anomaly corresponding, respectively, to the chiral ( "soft" ) and 
scale ("hard") symmetry-breaking contributions. Mixing depends on the baryonic sector to be described. 
For instance, in matter it is assumed that the glueball contribution decouples gradually with increasing 
density, reaching a critical point at which only the "soft" component is relevant. In such a case, the 
parameters can be fixed constraining the model by the bulk properties of finite nuclei and yielding the 
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values f x = 240 MeV and m x = 550 MeV. In section IV, we present the approximants to the solitonic 
solutions corresponding to the different sets of parameters reported elsewhere. Concerning the pion 
parameters, the accepted experimental values are f n = 93 MeV and m n = 140 MeV. 



III. PROPERTIES OF THE SINGLE SKYRMION WITH THE SCALAR FIELD 



In this section we analize the static solutions with nontrivial topology of the Euler-Lagrange equations 
of the model. We first present the exact asymptotic solutions, and then we propose suitable analytical 
continuations to these representations. 

Let us first derive the Euler-Lagrange equations from the Lagrangian Adopting the usual hedge- 
hog ansatz ^(| for the chiral field U = exp(ir • 71"//^) and assuming also spherical symmetry for the 
additional scalar field 9], that is 



U(r) = exp(ir • rF(r)), xir) = f x c ( r )> 



(5) 



we can write the mass of the single soliton in the Skyrme model as 



poo | f2 

E[F,C] = 4vr / r 2 dr 

Jo 2 



. f 2 

2 .;,. ) J 71 ( .2 



dF 



V 

dr J 



sin 2 F 



+ f 2 m 2 w C 3 {I- cos F) + J -+ 



f 2 
J _x 

2 



+ 2- 



dCY mi 



+ 



1 sin 2 F 
2^ 



=_j + _^ C >C-l/4) + i 



sin 2 F +2 (dFV 



+ 



(6) 



Variations of E[F,C] with respect to the profile function F(r) and the scalar field C(r) lead to the 
classical equations of motion. 

In terms of the dimensionless variable y, defined as y = ef n r, the coupled equations for the fields F 
and C read 



(y 2 C 2 + 2sm'F) 



d 2 F 



dy 



dC\ dF 



* +vC Hv)^ + 2sinFcosF {^l + 



dF 



2C 2 sin F cos F- 2 



sin F cos F I m 



and, 

d 2 C 2dC (f^ 2 
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dy 2 y dy V /; 



C 



dF\ 2 sin 2 F V 
+ 2 =- 



dy J 



y 



^) y 2 C 3 sinF = 0, 



> 2 / \ 2 



(7) 



■'^JfJ <? 2 (1 -cosF) - | rMnC'-O. 



(8) 



Finiteness of energy requires that U(r) tends to an arbitrary constant element of SU(2) at spatial 
infinity. Choosing U(r) — > 1 as r — > oo implies, for the chiral angle F(r), the boundary condition 



F(oo) = 0. 



0) 
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This condition compactifies the space of constant time surfaces, R 3 , to the three-sphere S 3 ; so the 
field U(r) defines a map from the compactified configuration space S 3 to the identity in the target space, 
SU (2), which is isomorphic to S 3 . The space of the distinct equivalence classes by homotopy corresponds 
to the group 

U 3 [SU(2)]=U 3 (S 3 ) = Z, (10) 

which means that the model has an infinite number of topological sectors, each one characterized by an 
integer- valued winding number. Z is defined as the topological charge q = J dx 3 Bo, which is related to 
the baryonic number current defined as 

B» = -L^ e ^Tr[(U^d u U)(U^d p U)(U^d a U)}. (11) 

As for the original Skyrme model, the equations of motion for F(y) and x(u) admit power-series 
solutions, of the form 

F{y) ~ K + F^ + Fsy 3 + 0(y 5 ), 

C(y) ~ C + C 2 y 2 + O(y 4 ), (12) 

which are valid for small values of y. The zero-order coefficient from the power-series representation 
of F(y) is choosed equal to ir in order to constrain the solution to the B = 1 sector. The coefficients 
F\ = —a and C$ cannot be determined by the perturbative analysis. 

By substitution of the power-series solutions 1)12(1 into Eqs.( |7| ) and (jSJ), the following recurrence rela- 
tions among the coefficients Fi, F3, Cq and C2 are found: 



C2 - C F? + - ^-2^^o ln Co, (13) 



2 KfxJ z 2 fl 6e2/2 

F3 = ~k 2 9 f ^ ( 2F yf* + 3m ' g o 3 + ^C,C 2 e 2 fl + 4C$F?e 2 fi). 

30 e z f4(C$ + 2FiY 

On the other hand, from the asymptotic behavior of the equations of motion it can be shown that 
the functions F and C reach their vacuum values at infinity following the functional forms: 

^ exp(-nWe/^ ^ ^ _ exp(-m x y/e/ 7r ) 

y y 

Our task in the next section will is to find approximate analytical representations for the functions 
F(r) and C(r), valid within the whole domain of definition of the radius r. 



Previous works 



IV. APPROXIMATE REPRESENTATIONS 

ing the equation of motion of the original Skyrme Model 
have shown the utility, for practical calculations, of implementing an approach that incorporates both 
asymptotic behaviors of the chiral angle F(r), near the origin and at large r, in one single representation. 
The rational approximants have proved to be suitable for this purpose. Moreover, it has been shown in 
Ref.[|| that starting from the well-known Pade approximant approach, a simple and systematic method 
to solve the Skyrme problem can be developed. In the following discussion we extend the analytical 
approach for the single Skyrmion with the scalar field, showing that the exponential behavior (|14j) can 
be easily incorporated into rational representations, providing satisfactory analytical continuations of the 
power-series representations We will begin by computing a sequence of Pade approximants using 

suitable ansatze for the functions F(r) and C(r) suggested from the previous study of the stucture of 
the exact solution. We continue the discussion by showing that the convergence to the exact solution is 
drastically improved when using a variational method to minimize E[F, C], supplemented by our ansatze 
for the fields. Then the attention is focused on analysing the convergence of the proposed solutions by 
use of some typical physical observables of Skyrmion-type models. 

For the sake of completness we briefly remind the reader that Pade approximants are rational functions 
used to give an analytic continuation of a power-series representation of a given function The Pade 
approximant P[m,n](u) °f order [M,N] to the series S(y) = J2n a nU n 1S defined as the ratio of two 
polynomials 

1 + Efc=i hv* 

where the free N + M + 1 parameters, otfc, are fixed so that the first M + N + 1 coefficients in the 
Taylor expansion of P[m,n] coincide with the series S(y) up to order M + N. 

In the present problem, modified approximants for the chiral angle F(r) may be built so as to match 
the exponential asymptotic behavior given in relation (|14j). 

Working with the dimensionless variable y defined previously, we may cast the approximate solutions 
for F in the form 

[ ) l + ZT=ihy k + b m+ iy m+1 eMm,y/ef n y { ' 

where the parameters and bk are to be fixed so as to reproduce the behavior of the exact solution near 
the origin. 

Note that the functional form that is proposed has the exact asymptotic behavior at infinity, that is, 
F(r) ~ exp(—m n r)/r. 
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The lowest order solution that can be built is 

7T 

1 + (a/7v)yexp(m n y/ef n )' 
The above approximant contains the asymptotic behavior, exp^—m^y/ef^/y, and agrees with the power- 
series solution near the origin up to order 0{y). It is possible to construct another approximant of the 
same order [0,1], but with an additional parameter b 2 , 

F (y) = — 7—^ , (18) 

1 + w V + b 2 yexp{m 7r y/ej n ) 
b 2 should be chosen so as to reproduce the exact asymptotic decay for large r. 

The higher-order approximants contain the information of additional terms of the exact expansion 
at the origin. Obviously, they improve the degree of approximation to the exact solution. For instance, 
keeping the ansatz (|16|) and fixing the coefficients by means of the Pade approximant method using 
equations (|12ll3j) . we find for the [1,2] approximant the following expression 

Mi, 2 ] = <l + aiy) (19) 

1 + hy + b 2 y 2 expim^y / f n e) 

with 

a nF 3 ttF 3 F 3 

ai = + — T - r , bi = — — -, 6 2 = i—r, (20) 

tt a z — iram-niej-K a z — iram n /ej~ n a — Tim^/ej^ 

and for the [2, 3] order approximant 

F r 2i 31 = tt(1 + a iy + a 2 y 2 ) 

1 + biy + b 2 y 2 + b 3 y 3 exp(m n y/ef 7T ) 

with 

ai = -- + 2(F| + aF 5 )(ir£ - a) /A, a 2 = [(a - (2a(F 3 2 + aF 5 ) - irF 2 ) - 2ir 2 F 3 F 5 ] /A, 

h = 2{Fi + aF 5 )(^-a)/A, b 2 = £F 3 « - 2a) - 2ttF 5 /A, b 3 = -2F 3 (F 2 + aF 5 )/A (22) 
where ^ = m^/ef^ and A = F3 [2-71-^3 + a£(iTt; — 2a)]. 

By construction, from the Taylor expansion of the Pade approximants ()19|) and 1)21(1 we recover the 
exact behavior near the origin [relations (|12|) ] up to orders 0(y 3 ) and 0(y 5 ), respectively. Note also that, 
as the proposed approximants satisfy the right boundary conditions, they yield to topological solutions 
of baryonic number equal to q = 1. That is, recalling Ea. (|lH) . 

dx 3 B (23) 



1 f°° dF sin 2 F 2 1 
-47ry ay 



2ir 2 J dy y 2 tt 



F{y) _ 



1. 
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The sequence of Pade approximants to the chiral angle F(y) given by equations (|17W21j) is shown in 
Fig.l together with the exact numerical solution. The approximate representations show convergence 
to the exact numerical solutions; unfortunately, we found that higher-order calculations in the Pade 
procedure lead to singular approximants being no longer relevant to represent F(r). 

In order to eliminate spurious representations with singularities in the real y axis one may proceed 
heuristically by adding restrictions to the numerator coefficients of the approximants. For instance, by 
simply setting = 1 for all k we will produce well-behaved and positive-definite approximants, as shown 
in Fig. 2. This procedure, which relies on simplicity, can be particularly useful in problems for which one 
seeks for an analytical description of Skyrmion-Skyrmion interactions such as in dibaryon models or other 
multi-Skyrmion models However the prescription of setting = 1 is too arbitrary and does not 
really select the solutions that minimize the energy functional E[F,C]. 

One can actually confirm this fact by studying the sequence of restricted rational approximants 
[m, m+1] for the chiral angle. The sequence of restricted Pade approximants is displayed in Fig. 2 where we 
can see a fast convergence to the [4, 5] order representation. As we should expect, the convergence to the 
numerical solution is very fast at the boundaries, however the representations underestimate in 1% to 2% 
in the 1 < y < 2 region. Higher order approximations slightly improve the agreement with the numerical 
solution but at the price of getting less attractive expression for using in phenomenological calculations 
from the model. Quality tests of the approximation are done further in this work by calculating physical 
quantities that frequently enter in the evaluation of the static properties of the nucleons. 

The low convergence to the numerical solution is a consequence of the fixing condition ai- = 1 which 
appears to be too restrictive. We will see further in this section that the convergence to the exact numeri- 
cal results is drastically improved and the arbitrariness in the approximate representations eliminated by 
taking the expansion coefficients in the ansatz (|16[) as free parameters and then minimizing by numerical 
methods the functional of energy E[F, C] given in Eq.ljBjl. 

To this end we first need to give a suitable ansatz for the scalar field C(y). One can notice that the 
asymptotic behavior of 1 — C{y) has the same functional form as the one corresponding to the profile 
function of the skyrmion so that analogous representations to those of F(y) should provide good results 
for the scalar field. Thus, we propose, 

Em k 

{ ' l + T!k=id k y k + d m+l y m +^Mm x y/eUy 1 ' 

Of course, we should expect a change on the dependence of the rational representations parameters 
on the power series coefficients as, in this case, the perturbative solution is given in even powers of y (see 
Eq-JUJl). 
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Our best representations of the exact solutions for F(y) and C(y) are obtained by keeping the func- 
tional forms (jlfij) . (|24j) and relaxing the conditions for the expansion coefficients. These are to be found 
by minimizing the functional of energy following a variational method. We give the first two approxi- 
mants in the sequence [m,m + 1] which are yet very good representations for the fields as it can be seen 
in Fig. 3. 



F[l,2] 
C[l,2] 
F[2,3] 
C[2,3] 



tt(1 - 0.097?/) 



(1 + 0.487y + 0.494y 2 exp(m 7r y/e/ 7r ) ' 
i (l- 0.729)(l + 2.003y) 

~ 1 + 2.511y + 0.409y 2 exp(m x y/e/ 7r ) ' 

tt(1 -0.342?/ + 0.103?/ 2 ) 

(1 + 0.353y + 0.171y 2 + 0.205y 3 exp( m7T y/ef n ) ' 
. (1 - 0.729)(1 -0.881y + 0.900y 2 ) 



(25) 



1 - 0.829y + 1.264y 2 + 0.095y 3 exp(m x y/e/ 7r ) ' 

We now show that the analytical representations are useful for phenomenological applications by 
making an explicit use of them and comparing the results with numerical calculations. We have used 
specific values for the set of parameters (m x , f x ) of the effective chiral Lagrangian, namely (720 MeV, 
240 MeV). From Ref.fl], we learn that these values are relevant for the discussion of nuclear matter. The 
analysis for other values of m x and f x , as those reported elsewhere Q], give analogous results. 

We have calculated the mass of the single soliton in the Skyrme model with the scalar field given by 
© as well as the moment of inertia of the soliton which, as it is well known, is found after performing the 
collective semi-classical expansion substituting U(r) by U(r,t) = A(t)U(r)A^ (t) in the Lagrangian, A(t) 
being an SU(2) matrix. For the model in question the collective transformation yields, after performing 
spatial integration, 



L = —E[F, C] + ATr^o^oA- 1 ] 



where E[F, C] is the soliton mass and the moment of inertia A is given by 



(26) 



A 



7T 



3jVe 3 



dy y 2 sin 2 F 



C 2 + 4 



dF 
dy 



sin 2 F 

y 2 



(27) 



For the mass of the soliton and the moment of inertia we found the following numerical values M = 1387 
MeV and A = 0.12//^, which we take as reference values when using the analytical representations for 
calculations. To further check the reliability of the rational representations we have also considered the 
axial-vector current which enters in calculations of matrix elements within the Skyrme model. Starting 
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TABLE I: Predictions from rational approximants obtained by minimizing Eq.JfJJl. Masses are given in units of 
/ir/e, and A in units of 1//^. 



[1,2] [2,3] [3,4] exact 


M 


71.10 70.86 70.86 70.87 


A 


0.110 0.116 0.116 0.116 


G 


0.67 0.63 0.64 0.66 



with Eq.Q one finds the axial- vector current to be 
( J a)£ = i f -f {£)^ [r^Utf -d^U) 



+ 



(28) 



16e 2 



TV 



Tr 



T a ,d v U*u] \daU ] U,d u U ] U 



so that the spatial integral of the axial current, expressed as a product of spatial and internal factors 
reads 



d 3 x(J A ) 



-1\ 



where the factor G is explicitly given by, 



a- 71 



dy y 



2 i dF . sin2F 



dy 



-GTi{T a AT j A 



sin2F f dF\ z sin 2 F dF sin 2 Fsin2F 
4 ( — ] +8 = — + 4 5 



(29) 



y 



dy J 



V dy 



(30) 



The values for the physical quantities obtained from the above representations are displayed in table 
I. We found a very good agreement between approximate and numerical solutions. 

We would like to stress that the phenomenological calculations, which involve a highly nontrivial 
nonlinear problem, get considerably simplified through the use of the proposed analytical approximate 
representations. Certainly, the most interesting feature of our representations is that they contain, in an 
explicit form, exact analytical properties of skyrmion solution for the Skyrme model with the scalar field 
problem as Eas. H12j) and (|14[) can be fully recovered from them. 



V. FINAL REMARKS 



The suggested approximants based on a Pade-like method provide satisfactory analytical represen- 
tations for the single soliton solution of an effective chiral Lagrangian, well-behaved under the scaling 
properties of QCD. The rational fractions proved to be well suited to incorporate the exponential behav- 
iors of the chiral and scalar fields. Therefore they can be used with reliability to implement an analytic 
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continuation of the exact power series solutions of the equations of motion near the origin. By this means 
we give an approach to a difficult problem of two coupled nonlinear differential equations where only 
numerical solutions are available. Besides, the analytical representations contain explicit information of 
the exact solution, which remains hidden in a numerical approach. Certainly, analytical representations 
for the fields simplify the analysis of physical situations described by the model and may be used to 
study the phase structure of the Skyrmion system modeling nuclear matter as well as in many other 
applications such as in the context of exotic baryons and monopole excitations where the chiral 
soliton models are the common background. 
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FIG. 1: Sequence of Pade approximants to the chiral angle F(y)/n and exact numerical solution for m x = 720 
MeV 
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FIG. 2: Sequence of constrained approximants to the chiral angle F(y)/ir and exact numerical solution for m x = 720 
MeV 
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FIG. 3: Sequence of approximants to F(y)/n and C(y) found by minimizing the functional of energy using rational 
representations for the fields. 



